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Neutron-TPC (nTPC) is a fast neutron spectrometer based on GEM-TPC (Gas Electron Multiplier~Time Pro- 
jection Chamber) and expected to be used in nuclear physics, nuclear reactor operation monitoring, and thermo- 
nuclear fusion plasma diagnostics. By measuring the recoiled proton energy and slopes of the proton tracks, 
the incident neutron energy can be deduced. It has higher n/y separation ability and higher detection efficiency 
than conventional neutron spectrometers. In this paper, neutron energy resolution of nTPC is studied using the 
analytical method. It is found that the neutron energy resolution is determined by 1) the proton energy resolu- 
tion (7£p/ Ep), and 2) standard deviation of slopes of the proton tracks caused by multiple Coulomb scattering 
(Cre(scateing)) and by the track fitting accuracy (74). Suggestions are made for optimizing energy resolution 
of nTPC. Proper choices of the cut parameters of reconstructed proton scattering angles (ew), the number of 
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ind the working gas help to improve the neutron energy resolution. 
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INTRODUCTION 


Neutron-TPC (nTPC) is a fast neutron spectrometer under 
research, and is expected to be used in the fields of nuclear 
physics, nuclear reactor operation monitoring, and thermo- 
nuclear fusion plasma diagnostics [1-3]. It is based on a 
GEM-TPC (Gas Electron Multiplier-Time Projection Cham- 
ber) using argon and hydrocarbon mixture as the working gas 
at 1 atm [4, 5]. A field-cage made of polyimide and copper is 
installed between the cathode and readout board to establish 
an effective volume, and a triple-GEM module works as the 
electron multiplier. A collimated neutron beam entering into 
the sensitive volume will scatter with protons of the working 
gas (Fig. 1). By measuring the energy deposited by the scat- 
tered protons (Ep) in the sensitive volume and the slopes of 
proton tracks (k), one can deduce the incident neutron energy 
by Eq. (1). 


Ey = Ep/ cos? 0 = Ep(1 + tan? 0) = Ep(1 +1/k°), (1) 


where Ø is the proton scattering angle, and k = tan(@ + 
7/2) = —1/tand. 


Fig. 1. Scheme of a neutron-TPC. 
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Thanks to TPC’s ability of 3D-track reconstruction and 
its large gaseous volume, neutron-TPC is advantageous in 
its higher n/y separation ability and detection efficiency 
than those of conventional neutron spectrometers, like the 
proton telescope system [2, 3, 6, 7]. By simulation based 
on the Geant4 software, neutron-TPC can reach a detection 
efficiency of ~0.1% and neutron energy resolution (FWHM) 
of better than 5% [1]. 

In this paper, we focus on deduction of analytical expres- 
sions of the neutron energy resolution of neutron-TPC, and 
discuss the optimization of neutron energy resolution based 
on simulations using Geant4 and Garfield [8-10]. Influences 
of different parameters on neutron energy resolution are stud- 
ied quantitatively, which helps optimizing the detector struc- 
ture and experimental parameters. Analytical expressions of 
the neutron-TPC’s energy resolution are deducted, and meth- 
ods to improve neutron energy resolution are suggested. 


Il. ANALYTICAL EXPRESSIONS 
A. nTPC energy resolution 


From Eq. (1), the neutron energy resolution (a, / En) can 
be derived by the resolutions of the proton energy (0 £p/ Ep) 
and the slope of proton track (a4 /k). 


(2) 


Assuming the incident neutrons are mono-energetic and only 
the protons with a certain scattering angle (0) are studied in 


040602-1 


202306.00250v1 


chinaXiv 


HUANG Meng et al. 


the data processing, then Ey is determined and the main task 
to improve the neutron energy resolution is to reduce standard 
deviations of the proton energy (a,) and the slope of the 
proton track (a). ). 


In the data processing, Ep is deduced by the number of 
electrons collected by the readout pads. As a result, the stan- 
dard deviation of the deduced proton energy (g,) is mainly 
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caused by the statistical fluctuation in the ionization process 
and the gain fluctuation of the GEM module. The slope k is 
deduced from the track reconstruction on the r-z plane. Two 
factors contribute to the slope uncertainty of the reconstructed 
proton track (a, ): one is the multiple Coulomb scattering of 
protons (cé(scatring)), another is the track fitting accuracy on 
the r-z plane (opci). Therefore, the neutron energy resolu- 
tion expression is the sum of three terms. 


or,\? _ (on, a —20 (scattering) \? q (2w a 
Ex] (E B+k B+k) ” 


In the following paragraphs, detailed analysis will be put 
on these three terms. 


B. Energy resolution of scattered protons 


First, assume there is no electron attachment in the work- 
ing gas and ignore the ADC (Analog-to-Digital Conversion) 
process. Then the proton energy deposited in the detector 
can be represented by electrons collected by the readout pads. 
Considering that the collected electrons come from two inde- 
pendent processes of the drift and avalanche, the number of 
electrons can be expressed as 

na=n-A, a) 
where n is the primary electrons ionized by the protons, and 
A is the amplification factor of the GEM module. Therefore, 
the proton energy resolution equals the relative standard de- 
viation of n4 


where F is the Fano factor, which is 0.2-0.5 for gases, W 
is the average ionization energy, and 03,/(A)? is the rela- 


Pcoutumb = [13.6 MeV /(8pe)] + 


where [p is the areal density of material, /3 is the velocity of 
the incident particle, p is the momentum of the particle, 2 is 


tive variance of A, which represents the amplification fea- 
ture of the GEM module. For the electron multipliers, the 
amplification factor obeys the Polya distribution [11], 


(a) Lema]: 
(6) 


where o is a free parameter determining the shape of distribu- 
tion. The relative variance of the Polya distribution is given 
by 03/(A)? = 1/(1 + a), and (a£p/ Ep)? can be given by 


(P -e (+e) (+r) 


From Eq. (5), one can find that, for a given neutron en- 
ergy, the proton energy resolution decreases with the slope 
of proton track, while the scattering angle increases and the 
scattered proton energy decreases. 


Utajo 


An PA= a 


(7) 


C. The o; (scattering) 


In measuring proton tracks, the unavoidable multiple 
Coulomb scattering affects accuracy of the reconstructed pro- 
ton scattering angle (9). At small angles, the multiple scat- 
tering angle £coutumb Obeys the Gaussian distribution, while 
at large angles, it obeys the Rutherford scattering. Here, 
€coutuny is defined as the angle between the initial direction 
of the recoiled proton and its final direction after a material 
of certain thickness. Based on the Highland Formula, one has 
the standard deviation of multiple scattering angles of inci- 
dent ions in material [12]: 


2(lo/Xo)!/?[1 + 0.038 In(lo/Xo)}, (8) 


the atomic number of the particle, and c is the light velocity. 
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The radiation length of material, Xo, can be calculated by the 
empirical formula [12]. 


Xo = 716.4A/[Z(Z + 1) n(287Z-"/)], (9) 


where A is the mass number of material, and Z is the atomic 
number of material. For a mixture or compound, the radiation 


length can be deduced by 
1 w 

fey, 10 

Xo xX o 


where X; is the radiation length of the i™ element and w; is 
the corresponding fraction by weight, 


Fig. 2. Schematic diagram of multiple Coulomb scattering. 


The standard deviation of proton recoil angles caused by 
multiple Coulomb scattering (Fig. 2) can be deduced [12]. 


o(scatering) = F ype / T 

= Tonin * T: 37? |e ay 

HBO pomm 
where z is the thickness of the material, ypiane is the lateral 
displacement of particle through the material in the projection 
plane. The slope of proton track and the scattering angle is 
related by k = —1/tan@ , sO o4¢scatering) Can be deduced 
from 074; scattering) 


iad 
(scattering) = 79(scatering) / Sin” O 


12) 
= 6 (scattering) (1 + k’). 
To prove correctness of the analytical expressions of 
(scatering)> a Monte Carlo simulation based on the Geant4 
code was carried out. In the simulation, 1000 protons events 
were run for every parameter setting, and their final track po- 
sitions were recorded to calculate (catering) The results are 
given in Table 1. The Geant4 simulation results are just a lit- 
tle smaller than those from the analytical expressions. There- 
fore, the analytical method is consistent to estimate the stan- 
dard deviation of the slope of the proton track caused by the 
multiple Coulomb scattering. 


RAT) 
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TABLE 1. Comparison of o9(scascing) from analytical expressions 
and Geant4 simulation 


Tas Thickness Ey Tan 
(cm) (MeV) (mrad) 
Analytical Simulation 
Ar 3 3 196 17.1£0.7 
5 11.8 9.8 £03 
10 3 28.8 
5 173 
Ar-C:Ho 5 3 15.6 
(50 : 50) 5 9.4 
10 3 29 
5 13.8 
D. The oa 


Another factor resulting in the slope uncertainty of the re- 
constructed proton track is the accuracy of the track fitting 
process, which can be expressed as oj.) in the neutron-TPC, 
the projections of proton tracks on the r-z plane are approx- 
imately linear due to absence of magnetic field. Therefore, 
the fitting function of proton track points (r;, 2;) is chosen as 
a linear function. 


eg bear, (13) 


For different proton scattering angles, the slopes of fitting 
functions c; are also different. For a certain proton scattering 
angle, one can get the corresponding uncertainties of param- 
eters co and c; based on the Least Squares Method [13]: 


Ve = (F"WyF)?, ag 


where V¢ is the covariance matrix of parameters co and c1. 


T o aaa 
Ve= | cov(criey) 02, (15) 
ln 
F=|: :J, (16) 
Tin 
and Wy is the weight matrix: 
oy 
o an 


Standard deviation of the fitting parameter c1, which is just 
«(ne)» Can be calculated by Eq. (18). 


2 _ 42 
Feet) = Fer 


poe 


(18) 
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where N is the number of points, r; and z; is the r- and z- 
coordinates of the fitting track point, and the oz, is the z- 
resolution of neutron-TPC at z = 2;. Obviously, apgr) has 
a close relationship with N and ø., that is, om) decreases 
with increasing N or improving o. 

Next, the most important step is to parse the z-resolution 
of neutron-TPC ø.. Generally, many factors will contribute 
to the z-resolution, including the intrinsic detector resolution, 
the gas properties, the tilt angle of the proton track with 2-axis 
etc. Equation (19) can be obtained to represent the influences 
of factors on the z-resolution. 


OF = Ofewctor + Faectonies + Talt angle + Faitusions (19) 


where oyeictor is the intrinsic detector resolution, which is re- 
lated to the readout pads layout, the GEM module setting, 
uniformity of the drift electric field, etc.; oetecuonies includes 
influences of electronic noise and analog-digital converting 
process: Ouin anple iS concerned with tilt angle of the proton 
track with z-axis, and can be expressed by Ref. [14]. 


Giittangle = d’ /(12 tan? 6 Net), (20) 


where d is the pad width, @ is the tilt angle (also the proton 
scattering angle), and Nur is the effective electron number: 
aittusion iS caused by the longitudinal diffusion of electrons, 
and can be expressed by Ref. [4]. 


Sutusion = DE2/ Nett, Qi) 


where D; is longitudinal diffusion coefficient of the working 
gas, 2 is the distance from the original position of ionized 
electrons to the readout board. 

Based on Eq. (6), the effective electron number Ny can be 
expressed as Eq. (22) [4]. 


p 1 = (5 + a) 
= = =a ; 
(1+ Gt) h Gata 

(22) 
where A is the amplification factor of the GEM module, N 
is the number of drift electrons per pad row, and a is the pa- 
rameter determining the shape of the Polya distribution. For 
N > 50, (1/N) is close to 1/(N), and Neg can be given 
approximately by Eq. (23) [15]. 


N) (j +a) . (23) 


Nese 
a 24a 


where (N) is the average number of drift electrons per pad 
row, and can be evaluated approximately based on the Bethe- 
Block formula [16]. 


(24) 


(25) 


—In(1 — (v/c)?) — (v/e)?], 
(26) 


where —(dE/dir)ion is the average energy loss per distance 
the protons traversing in the matter, d is the width of read- 
out pads, @ is the recoil angle of protons, W is the average 
ionization energy, 2 is the proton charge, N is the density of 
the matter atoms, Z is the atomic number of the material, v 
is the velocity of protons, £o is the vacuum permittivity, ¢ is 
the electron charge, mo is the electron rest mass, I = [pZ 
is the mean excitation potential (I) = 10 MeV), and c is the 
speed of light. Note that, the Bethe-Bloch formula is only 
suitable for the evaluation of the energy loss in material of 
charged particles with high energy (generally > 1 MeV for 
protons) [17, 18]. 


E. Evaluation of energy resolution of nTPC 


Based on equations above, the energy resolution of 
neutron-TPC under specific parameters can be estimated. If 
the incident neutron energy Ey is 5 MeV, the recoiled protons 
are produced 50 cm away from the readout board with a scat- 
tering angle of 30°, and the working gas is Ar-C)H, (50 : 50) 
at Latm, then the relative neutron energy resolution would be 
deduced as follows. E 

First, if œ = 0.5, F = 0.2, and W = 26eV (for 
Ar-C,Hg of ~ 50 : 50), then the value of (og/Ep)? = 
6.01 x 1076. Second, if the number of fitting points is 20% 
of the number of track points, then [—o4(catering)/(k® + 
A)? = 1.28 x 1074. Third, if the longitudinal diffusion 
coefficient D of electron is 231.4 um/cm?/? (based on the 
Garfield simulation for the drift field of 200 V/cm), the 
pad width is 2mm, and the z-resolution caused by detector 
and Geecwonies is 300m (based on experimental results), 
then —(dE/dzr)ion = 124MeV, Ner = 1.14 x 10%, and 
[okt / (X° + k)2 = 4.62 x 1075. As a result, the neu- 
tron energy resolution (o g, / En) = 1.3%, corresponding to a 
FWHM of 3.2%. It means that, if all the protons are of recoil 
angles less than 30°, the neutron energy resolution (FWHM) 
will be better than 3.2%. 


II. OPTIMIZATION OF NTPC ENERGY RESOLUTION 


For certain incident neutron energy and proton scattering 
angle used in the neutron energy calculation, the nTPC energy 
resolution is dominated by three parameters: the proton en- 
ergy resolution, standard deviations of t the proton track slope 
caused by multiple Coulomb scattering, and the accuracy of 
fitting procedure. In other words, the factors influencing 
the three parameters indirectly affect the reconstructed neu- 
tron energy resolution, such as the cut of reconstructed pro- 
ton scattering angle Gu, the number of fitting track points, 
and choice of the working gas. Therefore, optimization of 
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the nTPC energy resolution will focus on discussing these 
factors. Besides, the neutron detection efficiency should be 
taken into account, which is of importance, too, for a neutron 
spectrometer. 


A. Cut of reconstructed proton scattering angle (Yeu) 


As described in Section II A, for mono-energetic incident 
neutron beams, energy of the scattered protons increases, and 
the proton energy resolution improves, with decreasing pro- 
ton scattering angle. Besides, a higher proton energy leads 
to a longer proton track, which helps improve accuracy of 
the reconstructed proton scattering angle. Therefore, setting a 
proper cut of reconstructed proton scattering angle will prob- 
ably improve the neutron energy resolution. However, this 
will decrease the neutron detection efficiency. Considering 
distribution of the proton scattering angle cross-section (og) 
for fast neutron (< 10 MeV) [19]. 


og = (4oosin? @)/ 8 € [0,7/2], (21) 


where ap is the total cross-section of proton scattering. The 
neutron detection efficiency (7) can be derived as 


Bex hy 


n= m 


(28) 
4 [Oe _ sin(20eu) 
=m $. [fae Ew |, 


2 4 


where no is the intrinsic neutron detection efficiency of 
neutron-TPC. 

In order to demonstrate effects of the cutting angle 
method, a Monte Carlo simulation was carried out. The 
effective volume of neutron-TPC was based on a real TU- 
TPC (300mm x 500mm); the working gas was Ar-C,He 
(50 : 50); the drift electric field was set as 200 V/cm; the 
ring-shaped readout pad was of 2mm width; and the 5 MeV 
incident neutron beam was collimated along the z-axis. The 
number of fitting points (N) was 80% of the whole number 
of track points (Nn). The parameters of working gas, such as 
electron drift velocity (3.89 cm/s), and transverse and longi- 
tudinal diffusion coefficients of 293.3 and 231.4 um/cm?/?, 
respectively, were calculated by Garfield program. Proton re- 
coil and ionization processes were simulated by the Geant4 
code based on the physics list of QGSP_BIC_HP, and the en- 
ergy deposition of recoiled protons in the effective volume 
was recorded [20]. Then, electron drift, diffusion and track 
reconstruction were simulated by ROOT code using the fast 
MC method [21]. 

The results are shown in Fig. 3, one finds that the neu- 
tron energy resolution improves, but detection efficiency de- 
teriorates, with decreasing cutting angle (ex). With the pa- 
rameters set above, a neutron energy resolution (FWHM) of 
better than 5% and a detection efficiency of over 0.1% can be 
achieved at the cutting angle of 30°. Therefore, the essence of 
the cutting angle method is to compromise between the neu- 
tron energy resolution and the neutron detection efficiency. 


oar ows 
z P 
006 
H ~ oo È 
E ` i 
ja e m$ 
5 aa 
oo A f 
f m$ 
oor 
i 
oo omo 
C 0 2 0 50 
Cutting angle (degree) 
Fig. 3. (Color online) Neutron energy resolution and detection 


efficiency vs, cutting angle (zu). 


Note that, the simulation did not take into account the GEM 
module setting, the uniformity of drift field, electronics noise 
and the analog-to-digital conversion process, i.e. the parame- 
ters Fdetector ANd Getectronics in Eq. (19) were not included in the 
simulation. Therefore, this simulation just provides an ideal 
result approximately, and further studies concerning cdetsto 
and dzlectonies Will be carried out in the future. 


B. Number of fitting track points 


The number of fitting track points (N) affects standard de- 
viation of the slope of reconstructed proton track. On one 
hand, N determines the multiple Coulomb scattering angle of 
the reconstructed proton track. The multiple Coulomb scat- 
tering angle and the standard deviation of the slope increase 
with N. On the other hand, from Eq. (18), the increase of 
N decreases uncertainty caused by the track fitting process. 
Therefore, there exists an optimal value N, which leads to the 
smallest standard deviation of the slope, and the best neutron 
energy resolution. 
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Fig. 4. Neutron energy resolution vs. number of fitting track points, 
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Take the simulation with the parameters the same as in 
Sec. IITA, especially, feu = 30°. Considering the differences 
of proton track lengths at different recoil angles, a new vari- 
able R is defined as the ratio of the number of fitting points 
(N) and the whole number of track points (Nay). As shown 
in Fig. 4, at R = 0.2, the neutron energy resolution (FWHM) 
is the smallest (~2.2%). So, in data processing of nTPC, an 
optimized N is of great importance. 


C. Choice of the working gas 


In nTPC, the working gas plays both roles of ionizing 
medium and neutron detection medium. Generally, there are 
some basic principles of choosing a working gas of nTPC: 1) 
being of large enough detection efficiency for fast neutrons; 
2) ensuring proper operation of the GEM module; and 3) con- 
tributing to a good neutron energy resolution. This section 
will mainly discuss the influence of the working gas on the 
neutron energy resolution 

The characteristics of the working gas concerning the 
neutron energy resolution includes the transverse diffusion 
coefficient (Dr), the longitudinal diffusion coefficient (D4), 
the track length of recoiled proton (L), and the multiple 
Coulomb scattering angle of protons ({coulumb). First, the 
transverse and longitudinal diffusion coefficients affect 2- 
resolution of the nTPC, and small Dy and Di help achieving 
a good resolution of the slope of reconstructed proton track. 
Next, a longer L a larger whole number of track points (Na), 
which influences in turn choice of the number of fitting track 
points (N) and accuracy of the track fitting process. Finally, 
a smaller standard deviation of o eas leads to a smaller de- 
viation of the slope of the reconstructed proton track. As a 
whole, in order to reach a good neutron energy resolution, the 
working gas shall bear the following characteristics: small 
enough transverse and longitudinal diffusion coefficients, a 
proper proton track length, and a small enough standard de- 
viation of multiple Coulomb scattering angle of protons. 

Based on the Bragg-Kleeman rule, the track length in ma- 
terial can be calculated by Ref. [16] 


Li/Lo = mA; (1 Ay”), (29) 


where L is the track length, p is the density of the material, 
A is the atomic weight or effective atomic weight of the ma- 
terial, and the subscripts 1 and 0 denote the material studied 
and the air as reference, respectively. For compounds and 
complexes, the effective atomic weight (Aer) can be calcu- 
lated by AL? = D niA}/”, where n: is the number frac- 
tion of the i atom. Therefore, the density (p) and atomic 
weight (or effective atomic weight) of the working gas deter- 
mine the track length of protons in the drift chamber. So, a 
lower density and larger atomic weight (or effective atomic 
weight) leads to a longer proton track. 

From Eq. (8), the areal density (lo) and the atomic number 
(Z) of the material determine the multiple Coulomb scatter- 
ing angle for the same incident ion. Considering lo is product 
of density and thickness of the material, when the material 
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thickness is constant, the multiple scattering angle increase 
with the material density or the atomic number (Z). For the 
nTPC of constant N, the density (p ) and the atomic number 
(Z) of the working gas determine the standard deviation of 
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Fig. 5. (Color online) Neutron energy resolution as function of drift 
electric field, for the nTPC using different working gases. 


Argon and hydrocarbon mixtures are used as working gases 
for nTPCs, such as Ar-CH, (90 : 10), Ar-CH, (70 : 30), and 
Ar-C,H (50 : 50). Figure 5 shows the Geant4 simulation 
results of energy resolution of the nTPC using the three gas 
mixtures, as function of the drift electric field (Eain). In the 
simulation, the ratio of number of fitting track points and the 
whole number of track points was N/Nau 
cutting angle was feu = 30°. Note that the s 
not take into account the electric noise and ADC procedure, 
and the gas parameters in the different drift fields, such as 
electron drift velocity and diffusion coefficients, were cal- 
culated based on the Garfield simulation. From Fig. 5, at 
Eain = 300 V/cm, the neutron energy resolution of Ar-C)H, 
(50 : 50) reaches minimum for its relatively low transverse 
and longitudinal diffusion coefficient. With its relatively high 
hydrogen content, Ar-C,H, (50 : 50) is a good alternative as 
the working gas of nTPC. Studies will be done with more 
working gases. 


IV. CONCLUSION 


This paper mainly focuses on deduction of analytical ex- 
pressions of the neutron-TPC energy resolution. The neu- 
tron energy resolution (op, /E,) can be influenced by three 
factors: the proton energy resolution (o £p/ Ep), standard de- 
viation of the slope of proton track caused by the multiple 
Coulomb scattering angle ((scatering)) and by the fitting ac- 
curacy (c}(n:)). Based on the analytical expressions, one can 
estimate the neutron-TPC’s energy resolution under certain 
parameters, and optimize the resolution by proper choice of 
the cut of reconstructed proton scattering angle (eu), the 
number of fitting track points (N), and the working gas. The 
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setting and analysis of Aew and N are performed in the data 
processing, which have no direct relationship with the detec- 
tor’s structure, while the choice of the working gas is the key 
factor concerning optimization of the detector’s structure to 
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improve its energy resolution. From the Geant4 simulation 
results, Ar-C3H, (50 : 50) is a good choice as the working gas 
of nTPC, with higher detection efficiency for fast neutron and 
a better neutron energy resolution than those of other gases. 
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